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Abstract Adult males of the Papilio bianor complex were collected in Taiwan, 
Ishigaki island, the main island of Okinawa, and Hyogo. The samples represent 
members of certain subspecies of this species complex. Cell-free extracts of isolated 
testes and other arrayed tissues were analyzed by SDS-PAGE. After reading of the 
separated bands, the proportion of non-common bands out of the total was calculated 
and expressed by a dissimilarity index (DI). The average DI value within the 
populations was 0.017. On the other hand, that among localities was 0.153. The given 
DI values were 0.212 between Taiwan and Ishigaki, 0.148 between Taiwan and 
Okinawa, 0.128 between Taiwan and Hyogo, 0.160 between Ishigaki and Okinawa, 
0.116 between Ishigaki and Hyogo, and 0.155 between Okinawa and Hyogo. 

Papilio maackii is closely related to P. bianor complex, but it is a distinct species. 
When we take P. maackii as a control, the DI values between this species and each of 
the local populations of P. bianor complex were 0.145, 0.128, 0.102, and 0.178 for 
Taiwan, Ishigaki, Okinawa, and Hyogo, respectively. These values were almost equal 
to or rather smaller than the average value of 0.153 within the above localities of P. 
bianor complex. These results strongly support our previous findings that genetic 
divergence between local populations of P. bianor complex is on species level 
(KOMINAMI et al., 1987). Speciation is in progress in P. bianor complex, especially in 
island populations. 
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Introduction 


The swallow-tail butterfly, Papilio bianor, is distributed in some islands of Japan 
and Southeast Asia. Geographical variation has been noted with respect to the 
beautifully colored pattern on the wings, a pattern that varies from island to island. At 
present, P. bianor is divided into seven subspecies based mainly on morphological 
differences (FUJIOKA, 1975). Hybrid sterility, an indicator of divergence sufficient for 
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classification of an organism as an independent species, has not been confirmed 
sufficiently because of difficulties in breeding of the larvae. However, some results 
indicative of divergence within the P. bianor species have been obtained (AE, 1971 and 
1985 ; HAMA, 1976 ; WAKI, 1981). After introducing genetic analysis, ICHINOSE and 
HORIUCHI (1985) regarded an isolate on Ishigaki Island as a distinct species derived 
from the P. bianor complex. KAWAZOE and WAKABAYASHI (1976) also regarded that 
on the Okinawa islands as distinct. 

Our previous report (KOMINAMI et al., 1987) estimated the magnitude of genetic 
divergence among P. bianor subspecies using an index of genetic distance (D) (NEI, 
1972). The average D value was 0.967 among the subspecies. This is almost the same 
as the value of 1.009 between P. bianor and P. maackii, and rather higher than the 
value of 0.241 between Luehdorfia japonica and L. puziloi. In this report, further 
experimental results were obtained to confirm this genetic divergence among sub- 
species of the P. bianor complex. 


Materials and Methods 


Materials. Adult males of Papilio bianor CRAMER, 1777, were collected at the 
following four localities : (1) Nantou, Taiwan, in late July 1988. This sample repre- 
sents a stock of the subspecies P. bianor takasago NAKAHARA & ESAKI (SHIROZU, 1960). 
(2) Ishigaki Island, in early May 1988. The sample corresponds to a stock of the 
subspecies P. bianor junia JORDAN (KAWAZOE and WAKABAYASHI, 1976). (3) Nago, on 
the main island of Okinawa, in early May 1988. The sample corresponds to the 
subspecies P. bianor ryukyuensis FUJIOKA (FUJIOKA, 1975). However, KAWAZOE and 
WAKABAYASHI (1976) regarded stocks collected on the main island of Okinawa as a 
distinct species, P. okinawensis FRUHSTORFER, 1898. Thus, in this article the samples 
collected at Nago were tentatively treated as P. okinawensis (KOMINAMI et al., 1987). 
(4) Akazai, Hyogo Prefecture, in late May 1988. The sample corresponds to a stock of 
the subspecies of P. bianor dehaanii C. et R. FELDER (FUJIOKA, 1975). For each locality, 
20 males were randomly chosen and employed for description of protein variation. 

As a control experiment, Papilio maackii tutanus FENTON, 1881, a sibling species 
to P. bianor, were employed. Twenty adult males collected at Akazai, Hyogo 
Prefecture, in late May 1988 were used for this purpose. These butterflies were kept 
at — 80°C. Before use, the samples were dissected. The isolated testes and its arrayed 
tissues were homogenized in 1 ml of PBS buffer solution (pH 7.4) in a ceramic grider. 
After centrifugation at 15000 rpm for 5 minutes, the supernatant was taken and utilized 
as the sample solution. 

Sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). 
Twenty-five microliters of the sample solution was electrophoresed at a current of 15 
mA for 15 hours following the method of WEBER and OSBORN (1969). The spacer gel 
was 4% polyacrylamide ; and the separation gel, 12.5%. The running buffer was 45 
mM Tris-glycine (pH 8.3). The gel was stained with Coomassie brilliant blue R-250 or 
silver (Silver stain kit, Kanto Chem., Co.). 
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Fig. 1. Protein spectra obtained by SDS-PAGE of tissues from various subspecies of 
Papilio bianor complex : a. Papilio bianor takasago, collected in Taiwan ; b. P. 
b. junia, collected at Ishigaki; c. P. okinawensis (= P. b. ryukyuensis), collected 
at Nago ; d. P. b. dehaanii, collected at Akazai, Hyogo ; e. Papilio maackii, 
collected at Akazai, Hyogo. In each picture, the letters A to G above the lanes 
indicate individual males. Lane H corresponds to molecular weight standards 
whose weights are shown in daltons on the right margin. 


Two-dimensional gel electrophoresis. Seventy-five microliters of each sample 
solution was put into a glass syrinder containing a pH-gradient gel (pH 3.5 up to 9.5), 
and was electrophoresed at 400 V for 12 hours following O’FARRELL (1975). The gel 
was 3 % polyacrylamide containing 9.5 M urea and 2 % ampholine. After the run, the 





NII-Electronic Library Service 


The Lepidopterological Society of Japan 


4 Hirohiko Kominami, Megumi KAsAHARA, Kaoru SUMIYOSHI and Osamu YAMAGUCHI 


P. maackii 
AKAZAI 


0.014 


| 


P. maackii 


AKAZAI 


P. bianor P. bianor 
TAIWAN "0.020 AKAZAI 
P. bianor 
TAIWAN 





P. bianor PRLR P. okinawensis 
ISHIGAKI NAGO 
A 0.015 
P. bianor P. okinawensis. 
ISHIGAKI NAGO 


Fig. 2. Distribution of the dissimilarity values (DI) for all the possible combinations 
of the subspecies of P. bianor complex and P. maackii. 


column gel was overlayed on an SDS-PAGE gel and sealed with 1 % agarose. Then 
SDS-PAGE was performed. The final gel was stained by the silver stain method. 


Results 


1, Intrasubspecific variation in SDS-PAGE 

Twenty males were analyzed for protein variation for each subspecies. The mean 
number of distinguishable bands on SDS-PAGE gel was 44.3 per individual male of P. 
bianor takasago. Those numbers were 38.0, 48.3 and 45.0 for P. b. junia, P. okinawensis 
and P. b. dehaanii, respectively. That of P. maackii was 38.0. Using these records, a 
similarity index (SI) or dissimilarity index (DI) was calculated as follows (c.f. 

ODUM, 1971): 

SI = 2C/(A+B); DI =1 - SI 
Where, A and B are the number of bands detected in individuals A and B, respectively, 
and C is the sum of the bands commonly found in individuals A and B. This process 
was repeated for all the remaining males. Its average value was used for description 
of intrasubspecific variation (Dlintra). Dlinra values were 0.020, 0.019, 0.015 and 0.015 for 
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Fig. 3. Comparison of protein spectra obtained by SDS-PAGE of four subspecies of 
P. bianor complex and P. maackii. The two marginal lanes contained molecular 
weight standards whose weights are shown in daltons. 


P. b. takasago, P. b. junia, P. okinawensis and P. b. dehaanii, respectively. The value 
for the control P, maackii was 0.014. Thus, there is only a small variation, actually less 
than 0.02, within each subspecies. Protein spectra for each subspecies of P. bianor and 
for P. maackii are shown in Fig. 1. The results are also summarized in Fig. 2. 


2. Intersubspecific variation in SDS-PAGE 

Intrasubspecific variation was not so large as mantioned above ; so we examined 
intersubspecific variation (DIiter) next by comparing a couple of protein spectra 
randomly chosen. Calculation of SI and DI was made in the way as for intrasubspecific 
variation. The average Dlinter was 0.153. This value is ten times higher than that of 
the intrasubspecific difference (0.017). The given Dlinter values were 0.212 between P. 
b. takasago and P. b. junia, 0.148 between P. b. takasago and P. okinawensis, 0.128 
between P. b. takasago and P. b. dehaanii, 0.160 between P. b. junia and P. okinawensis, 
0.116 between P. b. junia and P. b. dehaanii, and 0.155 between P. okinawensis and P. 
b. dehaanii. All the possible combinations are represented in Fig. 2. On the other hand, 
the Dlinter values between P. maackii and each of the subspecies of P. bianor complex 
were 0.145, 0.128, 0.102 and 0.178 for P. b. takasago, P. b. junia, P. okinawensis and P. 
b. dehaanii, respectively, as shown in Figs. 2 and 3. The average value was 0.138, just 
slightly smaller than that between the subspecies themselves (0.153). 
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Fig. 4. Two-dimensional gel electrophoresis of the proteins from two subspecies of 
Papilio bianor complex. A. P. bianor takasago, collected in Taiwan; B. P. b. 
junia, collected at Ishigaki. Arrows indicate spots unique to a given subspecies. 


3. Unique pictures of two-dimensional gel electrophoresis for subspecies 

One individual male was randomly chosen from each of the subspecies, and its 
proteins were analyzed by two-dimensional gel electrophoresis (Fig. 4). Some unique 
spots were obtained, as indicated by the arrows in the figure. A supporting evidence 
for the protein divergence is observed. 


Discussion 


1. Correlation between Dliner value and genetic distance, D 

From the above experimental results, we can see the following tendency : any of 
the Dlinter values between the subspecies are equal to or not so different from any of 
the Dlinter values between P. maackii and a given subspecies of P. bianor complex. In 
other words, the degrees of genetic divergence between the subspecies are expected to 
be at the species level. 
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Fig. 5. Correlations between dissimilarity index of proteins and average hetero- 
zygosity of isozyme genes and between dissimilarity index of proteins and genetic 
distance of isozyme genes. The letter r indicates coefficient of correlation. 


KoMINAMI et al. (1987) quantitatively analyzed the degrees of genetic divergence 
found in P. bianor complex using the parameter of genetic distance, D. The D value 
between P. bianor and P. maackii was 1.009. This is quite near to the average value 
of 1.0 between species of organisms already examined (see the review of NEI, 1975). 
The average D value between local races or local populations is less than 0.01. The 
obtained D values were 0.869 between P. okinawensis and P. b. dehaanii, 0.737 between 
P. b. junia and P. b. dehaanii, and 1.268 between P. okinawensis and P. b. junia, 
respectively. The same tendency as noted above is thus seen with respect to genetic 
distance. Therefore, we can expect a positive correlation between these two indices. 
The results are shown in Fig. 5. The coefficient of correlation was 0.27, although there 
was no statistical significance (t = 0.561, P>0.05, d.f. = 4), due to the small number of 
the samples. There was also a tendency toward a positive correlation between Dlintra 
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values and the average heterozygosity at isozyme genes (KOMINAMI et al., 1987). 


2. Species diversity of P. bianor complex in island populations 

There is an extensive variation in the morphological characters in island popula- 
tions of P. bianor, e.g., color tone on wings (KAWAZOE and WAKABAYASHI, 1976), 
indices of wing traits (IZUMI et al., 1984), or the shape of male genitalia (ICHINOSE and 
HORIUCHI, 1985). Some of them are specific to certain island populations. When we 
make intersubspecific hybrids compulsorily, abnormality in sex ratio and hybrid 
sterility are usually observed (HAMA, 1976 ; WAKI, 1981). In addition to the above 
differences, the genetic divergence (KOMINAMI et al., 1988) and genetic dissimilarity 
(the present results) among subspecies are strongly suggestive of speciation in progress 
in P. bianor complex, especially in island populations. Any of the present subspecies 
of P. bianor complex may thus be reclassified into separate species. 
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H 要 


ヵ カラ ス ア ゲハ と その 近 紋 種 間 の タン パク レベ ル よ りう 解析 し た 分 化 
(小南 裕彦 ・ 笠 原 恵 ・ 住 吉 weUa 修 ) 


東南 アジ ア に 生息 する ヵ カラ ス ア ゲ へ 集団 中 に は 多量 の 遺伝 的 変異 が 保有 され て お り , これ ら は 特に 隔離 
され た 鳥 小 集 団 た に お いて は 遺伝 的 浮動 と 相まって 種 分 化 が 進ん で いる と こと が 期待 され る . この 問題 の 一 端 
を 調べ る た め に , OSE, 石垣 島 産 , 沖縄 本 島 産 , 本 州 産 の 各 亜 種 カ ラス アゲ へ 雄 20 匹 ず つ を 採集 し , 腹 
部 を 擦り 潰し た 全 上 澄み タン パク を SDS-PAGE お よび 二 次 元 電気 泳動 で 解析 し た . 合計 の バン ド 数 に 占め 
る 異な っ た ベンド の 割合 を 非 類似 度 指数 DD と し て 定量 化す る と , 亜 種 内 で は , 0.020, 0.019, 0.015, 
0.015 で あり , 対照 区 に と っ た ミヤ マカ ラス アゲ へ ハ 内 で は 0.014 で あっ た . 他方 , RBE LAB BER ER 
で は 0.212, 沖縄 本 島 産 , 本 州 産 と の 間 で は , 0.148, 0.128 で あっ た . BEBEL, MABE, 本 州 産 で 
は 0.160, 0.116 で あり , また 沖縄 本 島 産 と 本 州 産 で は 0.155 で あっ た . これ ら の 値 は , 別種 で ある ミヤ マ 
カラ ス ア ゲハ と の 間 の DI, すなわち 0.145, 0.128, 0.102, 0.178 と 殆ど 同 程度 で あり , これ ら の 亜 種 問 で 
分 化 が 種 の レベ ル ま で 進ん で いる 事 を 示し , 遺伝 的 距離 か ら 推定 し た 結果 と 良く 一 致し て いる . 
(Accepted November 7, 1990) 
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